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Abstract 

The effect of different stacking order of graphene multilayers on the electric field induced band gap 
is investigated. We considered a positively charged top and a negatively charged back gate in order 
to independently tune the band gap and the Fermi energy of three and four layer graphene systems. 
A tight-binding approach within a self-consistent Hartree approximation is used to calculate the 
induced charges on the different graphene layers. We found that the gap for trilayer graphene with 
the ABC stacking is much larger than the corresponding gap for the ABA trilayer. Also we predict 
that for four layers of graphene the energy gap strongly depends on the choice of stacking, and we 
found that the gap for the different types of stacking is much larger as compared to the case of 
Bernal stacking. Trigonal warping changes the size of the induced electronic gap by approximately 
30% for intermediate and large values of the induced electron density. 
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I. INTRODUCTION 



Graphene is a single layer of carbon atoms with hexagonal symmetry [1]. Multilayers of 
graphene can be stacked differently depending on the horizontal shift between consecutive 
graphene planes, leading to very different electronic properties |2J, e.g. to various band 
structures. 

A perpendicular electric field applied to bilayer graphene, with the AB stacking, can open 
an electronic gap between the valence and conduction bands |3l. This was shown indirectly 

n n 

by transport measurements [4], |5j. Later on, spectroscopic measurements confirmed the 



opening of a gap in the energy spectrum |6HlO|. The extension of these bilayer results to 
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in the case the perpendicular 



three and four layers of graphene was presented in Ref. 
electric field was realized by a single gate. It was found that such an electric field causes 
an energy gap which was found to be a nonmonotonic function of the gate voltage, and a 
reentrant opening and closing of the gap was predicted as a function of the electric field 
strength. In Ref. 12] the electronic band structure of the ABA-stacked trilayer graphene in 
the presence of back and top gates was invistigated. 

Recently, we generalized our previous results [ll| to the case when two, i.e. top and back, 
gates were applied to three as well as to four layers of graphene systems [13( . We found that 
due to the trigonal warping the obtained results do not exhibit electron-hole symmetry. A 
non-monotonic dependence of the true energy gap in trilayer graphene on the charge density 
on the gates was found. We also predicted an indirect gap with a non-monotonic dependence 
on the gate voltage. Four layers of graphene exhibit a larger energy gap as compared to 
the three layer system, which is a consequence of the fact that Dirac fermions are present 
in the AB stacked graphene multilayers in case of an odd number of layers, while for an 
even number of stacked graphene layers only charge carriers with a parabolic dispersion are 
present at low energies [14J . 

Using Raman spectroscopy measurements the graphitic flake thickness, i.e. the number 
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and 
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In Ref. 
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of graphene layers, can be obtained, as was demonstrated in Refs. 
tunable three-layer graphene single-electron transistor was experimentally realized showing 
a transport gap near the charge neutrality point. To our knowledge, up to now, no four layer 
system was studied experimentally. Electrical tunable energy gap systems are of interest 
from a fundamental point of view, but also for possible applications in electronics (e.g. for 
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transistors) and photonics (i.e. wavelength tuning of a laser). 

The electronic low-energy band structure of the ABC stacked multilayer graphene was 
studied within an effective mass approximation in Ref. Ill, with special attention to the 
Lifshitz transition, in which the Fermi circle breaks up into several pockets. 

In this paper we study the effect of different ways of stacking of multilayers of graphene 
on the electric field induced band gap by top and back gates. We limit ourselves to those 
stackings that have been found in graphite. The Bernal stacking (ABA), which has hexagonal 
symmetry, is common and stable, but some parts of graphite can also have rhombohedral 
one (the ABC stacking) 18| . The band structure of three and four layer graphene systems 
in the presence of a perpendicular electric field is obtained using a tight-binding approach, 
where we used a self-consistent Hartree approximation to calculate the induced charges on 
the different graphene layers. We found that the gap for trilayer graphene with the ABC 
stacking is much larger than the one for the ABA stacking, which was studied in Ref. [13 



Similarly for four layers of graphene the energy gap also strongly depends on the choice of 
stacking, and is smallest in case of Bernal stacking. When taking into account the circular 
asymmetry of the spectrum, which is a consequence of the trigonal warping, we found 
considerable changes in the size of the induced electronic gap for the considered systems at 
intermediate and high densities of total electrons induced on the layers. 

This paper is organized as follows. A short overview of our tight-binding approach with 
a description of the self-consistent calculation are given in Sec. HI] for the ABC stacked three 
layer graphene in the presence of top and back gates. The corresponding numerical results 
are also discussed here. In Sec. IHII we investigate four layer graphene with different stacking 
order in the presence of top and bottom gates. Sec. [IV] summarizes our conclusions. 



II. THREE LAYER GRAPHENE WITH THE ABC STACKING IN AN EXTER- 
NAL ELECTRIC FIELD 

We consider a system consisting of three layers of graphene with the ABC stacking, 
which is modeled as three coupled hexagonal lattices with inequivalent sites A, and Bi 
(i = 1, 2, 3 is the layer number) with A\ and A2, as well as A3 and B2 atoms on top of each 
other, as shown in Fig. [U We use the Slonczewski-Weiss-McClure (SWMcC) parameters, 
i.e. 70, 71, 72, 73, 74 of tight-binding couplings for bulk graphite. Within each layer the 
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interaction between nearest neighbor A{ and Bi atoms is described by the parameter 70. 
The strong coupling between nearest layers, i.e. between A\ — A2 and B2 — A3 atoms that 
lie directly above or below each other is given by 71, and the weaker nearest layer coupling 
by 73 (74), i.e. between sites B\ — B 2 and A 2 — B 3 (Bi — A 2 , Ai — B 2 , A 2 — A 3 and B 2 — B 3 ). 
The interaction between the next nearest layers (Bi — B 3 ) is determined by 72, as is shown 
in Fig. [H(a) and for comparison in Fig. QJb) we show the unit cell for the ABA trilayer. 
Using these parameters we compose the tight-binding Hamiltonian for three layer graphene 
with the ABC stacking, which has the form Ji| 
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Hn 


H13 
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D 3 
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(1) 



/ 



where the rows and columns are ordered according to atom A from layer 1, atom B from 
layer 1, atom A from layer 2, atom B from layer 2, etc, with the following two by two 
matrixes: 
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7o/* 



D 2 = D\ 
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32 
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where 



f(k x , k y ) = e lk * ao/V 3 + 2e- ik * ao/2 ' r3 coskya /2, 



(2a) 



(2b) 



(2c) 



(3) 



with a = 2. 46 A the in-plane lattice vector length. The Hamiltonian for the ABA stacking 



was discussed in Ref. 
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To control the density of electrons on the different graphene layers and independently 
the Fermi energy of the system, a top gate with a density of negative charges n t > (the 
electron excess density is positive) on it, and a back gate with a density of positive charges 
nt, < are applied to the trilayer (a schematic picture was presented in Fig. 1 of Ref. 1 13h . 
As a result a total excess density n = n\ + n 2 + n 3 is induced (n = n t + n&), with n\ 
the excess density on the closest layer to the top gate, n 3 on the closest layer to the back 
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gate, and n 2 is the excess density on the middle layer. In our model the top or back gate 
produces a uniform electric field E t f, = n t ^e/2eoK, and due to the induced charges on the 
graphene layers, in its turn, create fields Ei = nie/2e n with e the permittivity of vacuum 
and k the dielectric constant. There is a simple relation between the charge density on the 
gates and the voltage between the gate and the closest graphene layer: V t> b = en tt bd/2e n, 
where d is the distance from the gate to the closest graphene layer (usually d is equal to the 
oxide thickness, which is typically about 300nm). For our numerical calculations we use the 
value k = 2.3, which corresponds to graphene layers on SzC^. The difference between the 
charge densities induced on the individual layers of graphene creates asymmetries between 
the first and the second layers, as well as between the second and the third layers, which are 
determined by the corresponding change in the potential energies A lj2 and A 2;3 

Ai >2 (n) = a(n 2 + n 3 - \n b \), (4a) 

A 2 , 3 (n) = a(n 3 - |n 6 |), (4b) 

where a = 6 2 Cq/eqK, with cq = 3.35A the inter-layer distance. The Hamiltonian Eq. ([1]) in the 
presence of the top and back gates is modified and we have to add A li2 (n), and — A 2j3 (n) to 
the first and third layer on-site elements in Eq. ([1]). The tight binding Hamiltonian operates 
in the space of coefficients of the tight binding functions c( k ) = (c^, cb u ca 2 , cb 2 , ca 3 , cb 3 ), 
where = caX^) and cb 4 = cb^ k ) are the i-th layer coefficients for A and B type of 
atoms, respectively. The total eigenfunction of the system is then given by 

Ni N t 

*k(f) = J2 CA ^t^ +J2 c B^f{^, (5) 

i=l i=l 

with iV; the number of layers. By diagonalizing the Hamiltonian one can obtain the six 
coefficients (in Eq. (jSJ)) for fixed values of the layer asymmetries, from which we obtain the 
excess electronic densities on the individual layers: 

2 



m = - 

7T 



J ' dk x dk y (\c Al \ 2 +\c Bl \ 2 ). (6) 



The coefficients CA t and depend on the energetic band index. Here we are interested 
in the case when the Fermi energy is located in the band gap, and in order to find the 
redistribution of the electron density over the different layers in the valence bands one 
should integrate Eq. (]6]) over the Brillouin zone. The Fermi energy can be tuned into the 
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opened gap, when the magnitudes of the top and back gates are equal to each other but 
with opposite charges on them. The other case when the Fermi energy is located in the 
conduction or valence band was discussed in Ref. is] for the ABA stacked trilayer where we 
found that the obtained results do not exhibit electron-hole symmetry in the presence of 
trigonal warping. Using Eqs. (fi H4bj) and we evaluate the energy gap A at the i^-point 
and the true gap, A, self-consistently for a fixed total density n t + rib = rii + n 2 + n 3 (see 



Refs.a and 



131). 



In the following we will consider two cases. Firstly, we neglect all interactions except 
between the nearest neighbour atoms in the same layer and between the atoms of adjacent 
layers which are on top of each other, i.e. we put 72 = 73 = 74 = 7s = 0. This leads 
to a circular symmetric spectrum. In our calculations we used the parameter 70 = 3.12eV 
which leads to an in-plane velocity v = y/SjQd/Zh ~ 10 6 m/s, and for the interlayer coupling 
strength, we take 71 = 0.377 eV (see Ref. |20|), and for the interlayer distance cq = 3.35A 
Secondly, the full interaction case is studied where the interaction between the different 
atoms is expressed by the SWMcC parameters (72 = —0.0206, 73 = 0.29, 74 = 0.12, 75 = 
0.025), i.e. the effect of warping is included. 

Fig. [2] shows the band structure for trilayer graphene with the ABC stacking when charges 
on the top and back gates are opposite but equal in magnitude with — rib = n t = 10 13 cm~ 2 
when only 70,71 are taken into account (with k = 2.3), and the Fermi energy is located in 
the forbidden gap. Notice that there is conduction band - valence band symmetry around 
the Fermi energy, and the true gap A occurs away from the i^-point where the gap is 
A = 266meV > A = 195meV. For the ABA stacking for the case when only 70, 71 7^ the 
true gap is zero for all densities. 

When all the interactions between the different atoms are taken into account the surface 
of constant energy is no longer circular. In Fig. [3] we show the gap Ao at the .fT-point 
(dotted blue curve), and the true direct gap A (solid red curve) for trilayer graphene with 
the full interaction, as a function of the top gate density n t providing the back gate density 
— rib — n t- For comparison in the same figure we show also the corresponding results, A' 
(dashed red curve) and A' (dot-dashed blue curve) when only 70, 71 7^ 0. Notice, that for 
high densities {—rib = nt ~ 10 13 cm -2 ) the inclusion of the full interaction leads to a lowering 
of the true gap by 30%. It is interesting to note that similar values for the energy gap s and 



the relative difference between them was found for the case of bilayer AB graphene 



13j: the 
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true gap for the AB bilayer at — n& = n t ~ 10 13 cm~ 2 is 142meV when k = 2.3 and 198meV 
for the case of k — 1, when the full interaction is included. These results compare with 
169meV (re = 2.3) and 207meV (k = 1) for our ABC trilayer. 

This similarity becomes more remarkable, if we compare the layer densities induced by 
external gates for the ABA and ABC trilayers with the AB bilayer. For the ABA trilayer, 



when only a back gate was applied to the first layer 111 ], we found that ri\ = 6.1, n 2 = 3.2 
and ri3 — 1.2 at = 10 (in units 10 12 cm~ 2 ). The small amount of excess charges on the last 
layers was explained by the fact that the graphene layers screen the electric field and the 
layer asymmetries between the last layers, counted from the gate, are very small. The true 
gap for this system (A = 17meV) is smaller in comparison with the bilayer case, where for 
the latter n\ = 8.3 and n 2 = 2.8 (A = 97 .7meV) . Now, when only a back gate is applied to 
the ABC trilayer we find that the densities on the second and the third layers (counted from 
the back gate) are very close to each other: n 2 = n 3 ~ 2 at = 10 and ri\ = 6.24, which 
makes the ABC system distribution and the gap (with A = 117meV) similar to the AB 
bilayer ones. In Fig. [T^a) one can see that in the case of the ABC stacking there are never 3 
atoms stacked on top of each other, as in the case for the ABA. As a result the electric field 
(of the gate located near the first graphene layer for the ABC stacking) penetrates easier to 
the last layers inducing excess charges, while for the ABA stacking the electric field is much 
more strongly screened. 

When both gates are applied to the ABC trilayer graphene (when the full interaction is 
included) the excess charge densities at — rib = n t = 10, shown in Fig. HJ on the outer layers 
are —n\ = n 3 = 4.9 and in the middle layer is zero. Notice, that the excess charge densities 
on the bottom and the top layers are symmetric as in the case of the AB bilayer, as well 
as the gaps have also similar values. While for the ABA trilayer it was n\ = —3.84 and 
n 3 = 3.67 [ljij], and n 2 = 0.17 when —nj, = n t = 10; the inclusion of the full interaction 
in the ABA case makes the excess electron density in the middle layer different from zero, 
and it opens a small gap about 5meV. So, we see that the ABC system has a large gap, 
comparable with the AB bilayer one and behaves as a bilayer with shifted sheets, while the 
ABA opens up much smaller gap and is similar to the case of an AA stacked bilayer. 

In the case of the previous studied ABA trilayer [ll| the inclusion of trigonal warping 
leads to a non-monotonic behaviour of the gaps as a function of gate voltage, as well as 
a much stronger lowering of the true gap. Here, we found that the energy gaps for the 
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ABC stacked trilayer is much larger as compared to the case of the ABA trilayer. Fig. [5] 
shows a 3-D-plot and the corresponding contourplot of the highest valence band for three 
layer ABC stacked graphene near the .fT-point (A'-point is chosen as the origin, k = 2.3) 
for n t = —rib = 10 13 cm -2 . The lowest conduction band is again symmetric with the highest 
valence band just as in the case when only 70,71 7^ 0. Here, for the ABC stacking we fine 
three maxima, but did not find additional maxima as in the case of the ABA stacking jl3 | 
and SIS db result we do not observe an indirect gap. 



III. FOUR LAYER GRAPHENE SYSTEM IN AN EXTERNAL ELECTRIC FIELD 



Now, we consider the four layer graphene system, which can be arranged in many differ- 
ent ways as schematically shown in Figs. [6]^a-c). The tight-binding parameters 7$ and the 
interaction between the individual carbon atoms for all these cases are indicated in these 
figures. Four layer graphene is described by the Hamiltonian 

/ Di H12 #13 H u \ 
H21 D2 H23 H24 
H31 H32 D 3 H34 
\ H±i H42 #43 D A J 

where and Di with i = 1,2,3 are the matrix elements of the ABC trilayer given by 
Eqs. ( l2"aH2bj) and for the ABCA stacking we have 



H 



(7) 



H 



14 



H 



41 
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42 



H. 



24 




(8a) 



H 



13 



H. 



34 



Hi 



32- 



Da = Da 



while for the ABCC stacking these matrixes have the following form: 



(8b) 




75/2 





Da = D, 



(9a) 



(9b) 
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We consider a four layer graphene system with top and back gates, which induce a total 
excess density n = 111+112 + 113 + 714, where rij is the excess density on the zth layer as counted 
from the top gate. The corresponding change in the potential energy between consecutive 
layers is 

Ai (2 (n) = a(n 2 + n 3 + n 4 - \n b \), (10a) 
A 2j3 (n) = a(n 3 + n 4 - \n b \), (10b) 

A 3j4 (n) = a(n 4 - \n b \). (10c) 

By adding A 11 = A lj2 (n), A m = A li2 (n) + A 2>3 (r2) and A /y = A lj2 (n) + A 2 , 3 (n) + A 3>4 (n) 
to the on-site elements of the II, III and IV layer of the ABCA or the ABCC four layer 
Hamiltonian, respectively, we obtain the Hamiltonian in the presence of top and bottom 
gates. The eight coefficients cj^ = caX k ) and cb 4 = cbX k ), for fixed values of the layer 
asymmetries defined by Eqs. (llOaHlQcl) . can be obtained by diagonalizing the corresponding 
Hamiltonian. The electronic densities on the individual layers are given by Eq. (jHJ). The gaps 
Ao, A are evaluated self-consistently analogously as was done for the three layer system. 

The variation of the gap A at the if-point (dot-dashed red curve), the true direct 
gap A (solid red curve) and the true indirect gap (dotted blue curve) A kk > with the top 
gate density n t (n b = —n t ) for four layer graphene with the full interaction is shown in 
Fig. [7(a) for the ABCA stacked four layer graphene and in Fig. [7(b) for the ABCC stacking. 
One can see that for the ABCA stacking with full interaction the true direct gap is very 
close to the corresponding gap in the case of a trilayer with the ABC stacking, e.g., for 
n t = —n b = 10 13 cm~ 2 the true gap is about YllmeV for four layer graphene with the 
ABCA stacking and for the ABC trilayer it is 169meV. In Figs. 0(a) and 0(b) we present 
the layer densities for the ABCA and ABCC four layer graphene systems, respectively, and 
we include the curves for the densities in the ABC stacked trilayer for comparison in both 
figures (dashed curves). It is remarkable that the excess densities for the ABCA system on 
the outer as well as on the inner layers are symmetric. Notice, that the densities, shown 
in Fig. [U(a), on the outer layers for the ABCA are very close to the ABC trilayer graphene 
ones for all the values of n t : at — n b = n t = 10 for the ABCA n 4 = —n\ = 5.4, while for the 
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ABC trilayer graphene = —ri\ = 4.9 (in the units 10 12 cm~ 2 ). We see that also the ABCA 
four layer graphene behaves as the AB bilayer. The localization of the atoms (see Fig. [6]) 
can explain why the excess electron densities on the outer sheets for the ABCA system are 
larger than the densities for the ABC trilayer, and even in comparison with the AB bilayer 
densities (which has n 2 = —n\ = 3.7 at — rib = n t = 10). For the ABCA system there are 
never 4 atoms on top of each other, as is the case for the ABAB stacking. As a consequence, 
the electric field of the top gate, e.g. at A 3 (see Fig. Eta)) is screened only by the B 2 atom. 
Similarly, the field of the back gate at B 2 is screened only by A 3 atom. As a result, both 
these atoms feel the field of the top as well as the back gate, which leads to a decrease of 
excess charges on the inner layers (i.e. to a neutralization of these charges by the opposite 
gates). However, an outer layer (of the ABCA system) which is mainly charged by its closest 
gate does not feel the further located gate, since the latter is screened by the inner layers. 
In the AB bilayer the two sheets feel both gates, and consequently the excess charges (by 
absolute value) are less than in the outer layers of the ABCA. Due to this, also the gap 
for the AB bilayer is less (see the gap value in previous Section) than the ABCA one (for 
the same strength of the top and back gates). We see also, that the gap is large when the 
amount of excess charges in the inner layers is small (as it is for the ABCA system). We 
found also that for k = 1 the true gap is \83meV for the ABCA; the relative difference with 
the case of k = 2.3 is only about 5%. 

In contrast, for the ABCC four layer graphene the excess density on the third layer is 
larger, and on the fourth layer is smaller than the corresponding densities found in the case of 
the ABCA system. So, the increase in the excess densities as well as the density asymmetry 
in the inner sheets leads to a decrease of the gap. Also, the fact that the third and fourth 
sheets are not shifted, i.e. they have the AA stacking order, explain that in the ABCC four 
layer graphene the electric field opens up a smaller gap. In both cases we found a much 
larger gap (about 170meV for the ABCA stacking and 70meV for the ABCC stacking at 
n t = —rib = 10 13 cm -2 ) than in the case of the ABAB stacked four layer graphene [l3] (with 
bmeV for the same density). So, we see that from all the systems, considered in this paper 



and in Ref. 
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5, the Bernal stacking leads to the smallest gap. 
Figs. H^a) and (b) show 3-D-plots and corresponding contourplots of the highest valence 
and the lowest conduction bands near the fT-point (.ff-point is chosen as the origin, k = 2.3) 



in the case of rit = —rib = 10 13 cm 2 , for the ABCA and the ABCC stacking, respectively. 
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The conduction band for the ABCC stacking has a "Mexican hat" shape maxima and 
minima on a ring, as shown in the contourplot, e.g. there is a minimum at k x ao ~ —0.17 
and k y a = 0. In its turn the valence band has a local minimum between the two maxima 
at the plane k x a§ ~ —0.17. The asymmetry between the contourplots for the conduction 
and the valence bands for the ABCC (see Fig. 0(b)) leads to an indirect true gap. At low 
densities there is a true direct gap for the ABCC, but due to the overlap between the bands 
at different points in fc-space the indirect gap is negative as is shown in Fig. E^b), i.e. we 
have a semi-metal for low gate densities. For the ABCA we find only three minima in the 
conduction band and a symmetric valence band (see Fig. E](a)), analogously with the ABC 
trilayer case. For the ABCA system the indirect gap is smaller than the direct one at low 
densities, and they coincide at high densities. 



When finishing this paper we came aware of a recent preprint [2l| on the effect of an 
electric field on multilayers of graphene with different stacking. They used the simplest 
approximation where only 70,71 7^ 0. They argued that the inclusion of the other tight- 
binding parameters do not affect strongly the band structure and the true gap. However, 
our calculations show that the true gap can be changed by 30%. 

IV. CONCLUSIONS 

The effect of different stacking order on the electric field induced energy gap of three and 
four layers of graphene was investigated. For three- as well as for four-layer graphene the 
energy gap strongly depends on the choice of stacking, and we found that the gap is much 
larger than for the previously studied Bernal stacking. We found that the true gap for the 
ABC trilayer and the ABCA four layer graphene is comparable with the corresponding gap 
for bilayer graphene with Bernal stacking. The account of the circular asymmetry of the 
spectrum, which is a consequence of the trigonal warping, considerably changes the size of 
the induced electronic gap for the studied systems. 
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FIG. 2: (Colour online) The circular symmetric band structure of trilayer graphene with the ABC 
stacking order around the K-point when charges on the top and back gate are opposite but equal 
in magnitude, i.e. — n& = rij = 10 13 cm -2 , for the case when only 70,71 7^ 0. Horizontal dotted line 
is the Fermi level. The true gap A and the energy gap at the i^-point Aq are indicated. 
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FIG. 3: (Colour online) The dependence of the gap Ao (dotted blue curve) at the if -point, the true 
direct gap A (solid red curve) for the ABC trilayer graphene as a function of the top gate density 
nt providing the back gate density is — n& = rit. For comparison we show also the corresponding 
results, A' (dashed red curve) and A' (dot-dashed blue curve) when only 70,71 7^ 0. 
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FIG. 4: (Colour online) The charge density rtj on the different graphene layers for the the ABC 
trilayer with k = 2.3 and with the full interaction included, as a function of the charge density on 
the top gate n t with the back gate density n& = —fH- 
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FIG. 5: (Colour online) The highest valence band, with the corresponding contourplots for the 
ABC stacked trilayer graphene near the K-point (i^-point is chosen as the origin) with equal but 
opposite charges on the top and back gate when nt = —n^ = 10 13 cm -2 . The Fermi energy is 
located in the energy gap at E = 0. 




FIG. 6: (Colour online) Schematic of the couplings between the different (^4-white and B- black 
dots) sites for four layers of graphene for: (a) the ABCA, (b) the ABCC stacking, and (c) the 
ABAB Bernal stacking. 
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FIG. 7: (Colour online) The dependence of the gap Ao at the /f-point (dot-dashed red curve), the 
true direct gap A (solid red curve) and the true indirect gap (dotted blue curve) A^/ as a function 
of the top gate density nt for four layer graphene where we included the full interaction. The back 
gate density — n& = nt is the same (but opposite in sign) as the top gate. Results are shown for: 
a) the ABCA stacking, and b) the ABCC stacking. 




FIG. 8: (Colour online) The layer densities rii (solid curves) for the four layer system as a function 
of the charge density on the top gate rtj (providing — n\, = nt) when the full interaction is included: 
a) for the ABCA stacking and b) for the ABCC stacking. In both cases we added the results for 
the layer densities n\ (dashed curves) for the ABC stacked trilayer when the full interaction is 
included. 
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FIG. 9: (Colour online) 3D-plots and corresponding contourplots of the highest valence band 
(bottom figures) and the lowest conduction band (top figures) around the -fT-point (i'T-point is 
chosen as the origin k = 2.3) when nt = — nf, = I0 13 cm -2 for: (a) the ABCA and (b) the ABCC 
stacking of four layers of graphene. 
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